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Abstract 

We consider a tippe top modeled as an eccentric sphere, spinning on a 
horizontal table and subject to a sliding friction. Ignoring translational effects, 
we show that the system is reducible using a Routhian reduction technique. 
The reduced system is a two dimensional system of second order differential 
equations, that allows an elegant and compact way to retrieve the classification 
of tippe tops in six groups according to the existence and stability type of the 
steady states. 

1 Introduction 

The tippe top is a spinning top, consisting of a section of a sphere fitted with a 
short, cylindrical rod (the stem). One typically sets off the top by making it spin 
with the stem upwards, which we call, from now on the initial spin of the top. When 
the top is spun on a table, it will turn the stem down towards the table. When 
the stem touches the table, the top overturns and starts spinning on the stem. 
The overturning motion as we shall see is a transition from an unstable (relative) 
equilibrium to a stable one. Experimentally, it is known that such a transition only 
occurs when the spin speed exceeds a certain critical value. 

Let us make things more precise and first describe in some detail the model 
used throughout this paper. The tippe top is assumed to be a spherical rigid body 
with its center of mass e off the geometrical center of the sphere with radius R. 
In addition, we assume that this eccentric sphere has a mass distribution which is 
axially symmetric about the axis through the center of mass O and the geometrical 
center C of the sphere. The tippe top is subject to a holonomic constraint since 
we only consider motions of the tippe top on a horizontal plane. An important 
observation is that the system has an integral of motion, regardless of the model of 
the friction force that is used. This integral is called the Jellet J and is proportional 
to the initial spin no (equation (I2.7[0 . For the sake of completeness we mention 
here that a nonholonomic model (rolling without slipping) for the tippe top is not 
appropriate, since the equations of motion then allow two more integrals of motion 
(the energy and the Routh-integral) that prohibit the typical turning motion, cf. eg. 
[3 [5] . We will choose the friction force to be linearly dependent on the slip velocity 
of the contact point Q of the top and the horizontal plane. 

In earlier works [3 [15] conditions of stability for asymptotic states of the tippe 
top were retrieved by the Lyapunov function method starting with the Newton 
equations of motion. We take a different approach to the problem. Our goal is to 
describe the behaviour of the tippe top in terms of Lagrangian variables and thereby 
classify its asymptotic motions as function of non-dimensional (physical) parameters 



A/C, the inertia ratio and e/R, the eccentricity of the sphere and of the value the 
Jellet J. As main novelty, we show that by ignoring the translational velocity of the 
center of mass in the friction law, the Lagrangian formulation of the dynamics of the 
tippe top allows a restriction to a system on SO (3) which is amenable to Routhian 
reduction 2JJ. The reduced system allows a rather simple stability analysis of the 
relative equilibria based on the relationship of the value of the Jellet 's integral J 
and the tumbling angle 9, see below. Through the Routhian reduction we retrieve 
a simple stability criterium which leads to the same conclusions as in [7J [15] . 

The hypothesis of neglecting translational effects is similar to the one made 
by Bou-Rabee et al. [15]. They motivated it by reasoning that for all possible 
asymptotic states (the relative equilibria of the system) the velocity of the center 
of mass is zero [7], and, for this reason, in a neighborhood of these solutions the 
translational friction can be neglected. In |15) the dynamics of the spherical tippe 
top with small friction has been studied without such an approximation and a full 
analysis of the asymptotic long term dynamics of the system is given in terms of the 
Jellet and the eccentricity of the sphere and inertia ratio. Remarkable is the fact 
that the stability results obtained by means of the Routhian reduction procedure 
fully coincides with the results of [15 . This justifies a-posteriori the approximation 
assumption in the friction force and shows that it accurately describes the behavior 
of the tippe top. 

The main result of this paper is summarized in the following theorem, also 
compare with [TJ \W\ . 

Theorem 1 In the approximation of negligible translational effects, a spinning ec- 
centric sphere on an horizontal (perfectly hard) surface subject to a sliding friction 
is reducible with a Routhian reduction procedure The relative equilibria of the 

reduced system are precisely the steady states of the original system. They are purely 
rolling solutions and except for the trivial state of rest, they are of three types: 

(i) (non-inverted) vertically spinning top with center of mass straight below the 
geometric center; 

(ii) (inverted) vertically spinning top with center of mass straight above the geo- 
metric center; 

(Hi) intermediate spinning top, the top precesses about a vertical while spinning 
about its axle and rolling over the plane without gliding. 

The existence and stability type of these relative equilibria only depend on the inertia 
ratio ^ , the eccentricity of the sphere and the Jellet invariant J. In particular, 
six regimes are identified in terms of the Jellet invariant and only three exhibit the 
'tipping' behavior. 

It turns out that the vertical states always exist, and intermediate states may branch 
off from them. Qualitative bifurcation diagrams corresponding to the possible dif- 
ferent regimes are sketched in Fig. [T] 

The reader is referred to Sec. [3]- Sec. [4] for the details and specification of the 
parameter ranges. There are three main groups Group I, II and III determined by 
existence properties of intermediate equilibrium states (similar to the subdivision 
from ). Tippe tops of Group I may admit intermediate states for 9 > 9 C , where 
9 denotes the angle between the vertical and the symmetry axis of the top. Two 
subgroups are distinguished according to change in stability type of the intermediate 
states. Relevant is that for tippe tops belonging to this class the non-inverted 
position is always stable, so they never flip however large the initial spin when 
launched under an angle 9 close to 0. Tippe tops of Group II may admit intermediate 
states for all 9, and they show complete inversion when the initial spin is large 
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Figure 1: Bifurcation diagrams of relative equilibria in function of the Jellet invariant J. Solid 
black branches correspond to stable relative equilibria, while dashed black branches correspond to 
unstable ones. The vertical states 8 = and 6 = n always exist, from which intermediate states 
(with < 9 < it) may branch off. 



enough. Tops of Group III tend to flip over up to a certain angle 9 C < n when 
spun rapidly enough. Since stability results are often in the literature expressed 
using the 'initial' spin no of the tippe top, one can read the J 2 in the figures as 
Uq. We anticipate to further results by noting that the instability inequalities are 
independent of the friction coefficient unless it is zero. 

The structure of the paper is as follows. In Section [2] our model is described and 
the equation of motion are given according to the Lagrangian formalism. After 
introducing the Jellet integral of motion, the Routhian reduction is performed in 
Section 12.21 The steady states of the system are then calculated and their sta- 
bility type is determined, yielding a tippe top classification in six groups which is 
summarized in Section [4. II 



2 Equations of motion 

As is mentioned in the introduction, we consider the eccentric sphere model of 
such a top, see Fig. [2] That is, we consider a sphere with radius R whose mass 
distribution is axially symmetric but not spherically symmetric, so that the center 
of mass and the geometric center do not coincide. The line joining the center of 
mass and the geometrical center is an axis of inertial symmetry, that is, in the plane 
perpendicular to this axis the inertia tensor of the sphere has two equal principal 
moments of inertia A — B. The inertia moment along the axis of symmetry is 
denoted by C and the total mass of the sphere is m. 

The eccentricity e is the distance between the center of mass O and the geometric 
center C of the sphere, with < e < R. The point Q is the point of contact with the 
plane of support. We assume that an inertial (laboratory) frame Mxyz is chosen, 
where M is some point on the table and the z-axis is the vertical. Let us denote 
the unit vectors along the axis of the reference frame Ox, Oy, Oz fixed to the body 
by respectively e x ,e y ,e z . The coordinates of the center of mass arc denoted by 
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z + n 




Figure 2: Eccentric sphere version of the tippe top. R is the radius of the sphere, the center 
of mass O is off center by e. The top spins on a horizontal table with point of contact Q. The 
axis of symmetry is Oz and the vertical axis is Oz, they define a plane El (containing OQ) which 
precesses about Oz with angular velocity ip. The height of O above the table is h(8). 



r o = ( X >y> z )-Mxyz- A second reference frame is denoted by Oxyz, and is defined 
in such a way that its third axis is precisely the symmetry axis of the top and the 
y-axis is perpendicular to the plane II through the z- and z-axes (see Fig. [5]). Again 
we denote the unit vectors along the axis of this reference frame by e x ,e y , eJ3- 

Let (9, <p, ip) be the Euler angles of the body with respect to the inertial frame, 
Fig. [2j chosen in such a way that (i) the vertical plane II is inclined at tp to the 
fixed vertical plane xz and precesses with angular velocity ip around the vertical Oz; 
(ii) the angle 9 is the angle between the vertical Oz and the axle Oz of the top; 9 
causes the nodding (nutation) of the axle in the vertical plane II; and (iii) the angle 
ip orients the body with respect to the fixed-body frame, ip is the spin about the 
axle. 

As it was pointed out before, the tippe top is constrained to move on a horizontal 
plane. This holonomic constraint is expressed by z = R — ecos9 = h(9). We 
assume throughout the paper that the only forces acting on the sphere are gravity 
G = —mge z and a friction force F exerted at the point of contact Q of the sphere 
with the plane. It is now immediate to write down the Lagrangian for the tippe 
top: 

L = \ (m(x 2 + y 2 ) + (e 2 m sin 2 9 + A)9 2 + A sin 2 9ip 2 + C(ip + <p cos 6>) 2 ) 

- mg(R - ecosd), (2.1) 

where g is earth acceleration. 

This function is defined on the tangent space of the configuration manifold 
M = IR 2 x SO(3). In order to obtain the equations of motion, it only remains 
to define a suitable friction force F. The Lagrangian equations of motion for the 
tippe top then read: 

<-> 

where q l represents one of the coordinates (x, y, <p, 9, ip), and Q F — Qf dq l is a one- 
form on M representing the generalized force moment of the friction force at the 
point of contact. It is defined by, with F = Ti,f-\-R n e z the orthogonal decomposition 
of F: 

Q F = R/ ■ e x dx + R/ ■ e y dy + (q x R/) • (e y d9 + e z d(p + e z dip). 



1 ln |16| the origin of the reference system attached to the body is in the center of the sphere 
C, and not in the center of mass O. 
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Modeling the friction force One typically models the friction force F = + 
R n e z to be proportional to the slip velocity of the point of contact vq. We denote 
by R„ = R n e z the normal reaction of the floor at Q, which is of order mg, and 
R/ = Fxe x + Fye y is the (sliding) friction which opposes the slipping motion of the 
body. The fact that the sliding friction opposes the slipping motion is expressed by 
R/ • v q < 0. We adopt a viscous friction law [TJ [5J |H1 HO] and assume that 

R/ = — /ii? n vg. (2.3) 

Here n is a coefficient of friction with the dimension of (velocity) -1 . It now takes 
a few tedious computations to arrive to the coordinate expressions for the force 
moments of the friction force. The coordinates of the point of contact OQ := q 
are Q = (xq, 0, ZQ)oxyz = (i?sin#,0,e — Rcos6)o X yz- The velocity of the point of 
contact Q equals 

v Q = v + wxq ) (2.4) 

vq = (x, y, h 1 (6)6) Mxyz is the velocity of the center of mass. A coordinate expression 
for the angular velocity is given by 

u = — p s'm(6)e x + 6e y + ne z , where n := tp + pcos(6) (2-5) 

n is the spin (that is, the component of u> about Oz). The generalized force moments 
now read: 

Q x = — /ii?„(x — sin ip 6(R — ecosff) + cos tp sin 6 (Rip + tip)) 

Q y = — fiR n (y + cos p 6 (R— t cos 6) + sin p> sin 6 (Rip + tp)) 

Qe = — nR n (R — ecos6')(cos py — sinpx + (R — t cos 6)6) 

Q v — — iiR n t sin 6 (cos px + sint^y + sin6(tp + Rip)) 

Qtp = — /J,R n R sin 6 (cos px + sin py + sin 6 (Rip + tp)). 

For the sake of completeness we write an explicit expression for the normal compo- 
nent of the reaction force R n , which can be determined from Newton's law for the 
center of mass of the sphere: 

R n (6, p, 9, p, iP, x, y) = 9 + §2k " + ^ COSi9) - + ^ COS 9),A) (2.6) 

1/m + h'/A[— h/i(sinpx — cosyy — 6h) + h'] 

To conclude this section, we briefly discuss other models for the tippe top. The 
eccentric sphere model does not accurately model the contact effects of the tippe 
top stem. However, it does describe the fundamental phenomenon of the over- 
turning. For the sake of completeness we remark that as soon as the top rises to 
spin on its stem, it behaves as a 'normal' spinning top with rounded peg, we refer 
to [HI Chapter 6] for a satisfactory introduction to this topic. Note that our model 
also does not describe the peculiar 'Hycaro' tippe tops by Prof. T. Tokieda [15] , 
which need a non-axisymmetric asymmetric mass distribution. These tippe tops 
have a 'preferred direction' meaning that the top would flip over only if spun in the 
preferred direction with a certain initial spin, and, no matter what the initial spin 
is, it would just continue rotating around the rest position when spun the other 
way round. We refer to [4 for a detailed analysis when elasticity properties of the 
horizontal surface and tippe top are taken into account. Their model allows for 
jumps of the tippe top on the horizontal surface (that we assumed to be rigid). 
Finally, we mention that we do not debate over the issue of whether transitions 
sliding-rolling and rolling-sliding occur in the motion of an eccentric sphere on a 
flat surface. We chose to concentrate on the sliding model only, because we were 
interested in capturing the 'overturning' phenomenon which cannot occur under the 
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non-holonomic constraint of pure rolling. We refer the interested reader to [TJ [5] 
for a discussion of the topic and to [2 [13l [19] for an analysis of the motion of the 
rolling eccentric sphere also called the Routh's sphere. 

2.1 Constants of motion: the Jellet invariant 

It was first shown by Jellet [B] by an approximate argument, and later proved by 
Routh [17] that the system, even if dissipative, has a conserved quantity: 

J = — L ■ q = const, (2.7) 

where L is the angular momentum of the tippe top about the center of mass. 
We prove this by using Euler equations which govern the evolution of the angular 
momentum L = qx F. The total time derivative of J then becomes: 

J = L • q — L • q 
= - (Auj -{A- C)(uj ■ e z )e z ) ■ (ew x e 2 ) = . 

Straightforward calculations show that Jellet 's constant can be written as 

J = Cn(Rcos(6) - e) + AtpRsm 2 (6) . (2.8) 

We emphasize once more that the Jellet 's constant is an exact constant of motion 
for the tippe top whether or not there is slipping and independent of the expression 
for F. As we will explain later, it is this constant that to some extent controls the 
motion of the spinning top. Indeed, it allows a Routhian reduction procedure (see 
Sec. 12. 2[ ). resulting in relatively simple reduced equations from which we are able 
to recover in full detail the stability properties of the steady states. In the specific 
case that 8 = 0, the Jellet is proportional to hq = n\g = o, the spin about the z-axis. 
Since one typically sets off the tippe top at an angle 6 ~ 0, one can say that the 
Jellet is proportional to the initial spin no. Note that the spin at 9 = tt has an 
opposite sign to the initial spin no, meaning that, relative to a body fixed frame, 
the spin is reversed when the tippe top fulfills a complete inversion. 

There is a rotational symmetry for which the Jellet is the associated first integral. 
The action of S 1 on R 2 x SO (3) can be defined as a simultaneous rotation about 
ez over the angle i?£ and about k over the angle — e£, where £, £ S 1 (see also [H]). 
Noether's theorem is applicable in this situation since the work of the friction force 
at the point of contact vanishes under this action. 

Note that the total energy of the spinning top is E = T + V is in general not 
conserved. Here T is the kinetic energy with its rotational and translational part, 
V = mgh{6) is the potential energy. The orbital derivative of E is 

j f E = v Q -K f <0, (2.9) 

which is negative semi-definite and vanishes if and only if vq vanishes. Observe that 
E(t) decreases monotonically and hence is a suitable Lyapunov function^, see [7]- 
From (|2.9|) it follows that dissipation is due to friction. 

2 E(t) is analytical, therefore it is either strictly monotone or a constant. The energy E is 
constant only if vq = 0. Note that E being Lyapunov [7] implies that the limiting solutions for 
t — > oo are solutions of constant energy. 
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2.2 Routhian reduction 



It turns out that, if we consider an approximation of the friction law, the resulting 
generalized force moments assume a form that allows us to apply a Routhian re- 
duction procedure, see |22| and appendix [X] In turn, using the reduced equations 
we are able to study in full detail the stability properties of the tippe top which 
confirm the results obtained in [15] , and also recover those of [16j [20] . 

We now ignore translational effects in the friction force, i.e. we assume that all 
terms in Qg, Q v , 0i/> containing x and y are neglected. Typically this approximation 
is justified by noting that for all steady states the velocity of the center of mass is 
zero, and that in a neighborhood of the steady states it can be neglected. In our 
situation, it allows to restrict ourselves to a system on SO (3) which is reducible 
using Routh's procedure. It is easily seen that within this approximation, if we 
study the Lagrangian system on 5*0(3) determined by 

V = \ {{e 2 msm 2 9 + A)9 2 + Asm 2 9tp 2 + C {ip + ip cos 9) 2 ^) 

—mg(R — ecos#), 
Q! = -fiR' n (R- e cos 9) 2 9d9 - fiR' n e sin 2 9{ep + Rip)dp 
-fxR^Rsin 2 6(Rip + eip)dip, 

with R' n (8, 8, ip, ip) = R n (9, <p, 9, ip, ip , x = 0, y = 0), then we essentially study the en- 
tire approximated system. Indeed, any solution (tp(t),9(t),ip(t)) to this Lagrangian 
system will determine the remaining unknowns (x(t),y(t)) as solutions to the fol- 
lowing system of time-dependent second order differential equations: 

mx=Q x (x,9{t),cp(t),iP(t)), my=Q y (y,6(t),tp(t),ip(t)). (2.10) 

Our next step in the reduction procedure is to consider the Lagrangian system L' 
on 5*0(3) with generalized force form Q' and perform a simple coordinate transfor- 
mation, determined by 

(9, ip, ip) H> (9, Tp = ep + Rip, c = Rtp — eip). 

The Lagrangian L' and the force Q' then become 

L' = \ (V 2 m sin 2 9 + A)9 2 + ^f^, (eip + Rc) 2 + 

( e 2^.2)2 ((R + e cos 9)Tp + (Rcos9 — e)c) 2 ^ — mg(R — ecos6*) 
Q' = -fj,R' n ((R - e cos 9) 2 9d9 + sin 2 9ipdTp) , 



where it is understood that R' n is a function of 9, 9, Tp and c. The main reason 
for writing V and Q' in this form is the fact that c is a cyclic coordinate and 
that Q' c — 0. Indeed, recall that Jellet integral was associated to the symmetry 
determined by the vector field Rd v — ed^, or in the above introduced coordinate 
system by the vector field d c . In particular, we have made the symmetry generator 
into a coordinate vector field, and this ensures that we can apply the Routhian 
reduction procedure ([22] and appendix E]) , provided the coefficients of Q 1 do not 
depend on c. This is the case since we neglected the terms in the velocity of the 
center of mass. The conserved quantity associated with the cyclic coordinate is 
precisely the Jellet integral: 

= J/(e 2 + R 2 ). 
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The latter equality implies that 

. _ J(e 2 + R 2 ) - (RAe sin 2 9 + C(R cos 9 - e){R + e cos 9))Tp 

R 2 A sin 2 9 + C(R cos - e) 2 ' ' ^ 

The Routhian reduction procedure defines a Lagrangian system 

TZ = L' — Jc/(e 2 + R 2 ) 

with two degrees of freedom (9,Tp) (here every instance of c in 1Z is replaced us- 
ing (|2.1ip ). The reduced equations of motion are then given by 

!(!|) -it - <*--"*-'* <«3) 

where it is understood that (|2.11l) is used to eliminate c in R' n . It takes rather 
tedious but straightforward computations to show that 1Z can be written as 



K = T 2 + T 1 -W, 



where 



T 2 = Tgg{9)9 2 +T^{9)Tp 2 



l , 2 .2/. «^2 1 AC sin 2 9 ^2 

= -(em sin + A)9 H = TB, 

2 V ' 2(R2Asm 2 9 + C(Rcos9-e) 2 ) 

rp _ J (RcAsin 2 0+C(R+ecos6)(Rcos6-e)) — 

11 ~ (R 2 +e 2 ) R 2 Asin 2 9+C(RcosO-t) 2 ^' 

w = l mAsi^e+c{R COS e-er ~ m 9^cos9. 
The function W is also called the effective potential. 

Remark 1 The above defined Routhian function 1Z is not globally defined on the 
sphere. In order to provide a globally defined system of differential equations for the 
reduced system, we need to extract the term T\ from the Routhian and consider it 
as a gyroscopic force (see e.g. \ '2 Vp. 

Remark 2 Observe that the effective potential W , obtained through reduction, co- 
incides with the effective energy on a Jellet's level surface as it has been used in 
Hf, Sec. 3. 



Note that TZ and Q' are independent of Jp. This residual symmetry does not lead 
to a conserved quantity (the friction does not vanish for dip). This symmetry is 
due to the approximation we carried out in the previous section; it is not present 
in the original system (L, Q) or (|2.2[) . It leads however to a zero eigenvalue of the 
linearized system at equilibrium points, see Sectional 

3 Steady states 

The equilibria of the reduced Routhian system (12 . 12[) and (|2.13p are determined by 

= 0,^=0 and 8W/d9 = 0, (3.1) 
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i.e. they are the critical points of the effective potential (note that if 9 = tp = the 
components of the force vanish). Equation (|3.ip is satisfied if (i) sin 9 — or, if (ii) 



f(j\™ S e) = _4 5 -((^-l)cos0+^)~(^sin 2 




(3.2) 



Solutions to (i) are 9 — 0, n and give the vertical spinning states. Solutions to (ii) 
only occur if 



If this condition is satisfied, solutions to (ii) are the so-called intermediate states. 
The existence condition only depends on A/C and e/R, and will determine in our 
classification the three main groups I, II and III, as it was proposed in [20] , The 
values for 9 for which (ii) is satisfied will depend on the Jellet, A/C and e/R. 

The vertical spinning states correspond to the periodic motion of the tippe top 
spinning about its axle (which is in vertical position) either in the non-inverted 
position (9 = 0) or inverted position (9 = n). The intermediate states correspond 
to those relative equilibria in which the tippe top shows in general quasi-periodic 
motion precessing about a vertical while spinning about its inclined axle rolling over 
the plane without gliding (observe that the intermediate states correspond to the 
tumbling solution of [7])- 

The condition (|3.3[) for existence of intermediate states leads to a first classifi- 
cation of tippe tops into three groups. 

- Group I [{A/C — I ) < —e/R}: the tippe tops belonging to this group do not ad- 
mit intermediate states in an interval of the form [0,6> c [, where 9 C is determined by 
(A/C - I) cos^c + e/R = 0. 

- Group II [—e/R < (A/C— 1) < e/R}: intermediate states may exist for all 9 G]0,7r[. 

- Group III [(A/C — 1) > e/R]: tippe tops belonging to this group do not admit 
intermediate state in an interval of the form }9 C , tt[, where cos^c = (e/R)/(l — A/C). 

In the following section we refine this first classification taking into account the 
stability type of the steady states and their bifurcations, explaining and giving the 
details of the J 2 versus 9 diagrams in Fig. Q] from the introduction. 

We anticipate that the subdivision in subgroups according to a change in sta- 
bility type of the intermediate states is based on the simple observation that they 
lie on a curve f(J 2 , cos 9) = in the (J 2 , #)-plane, and, denoting by 82 f the partial 
derivative to the second argument of /, a bifurcation point for intermediate states 
is characterized as the point where 



Note that the relation c?2/( J 2 , cos 9) = is essentially the same basic relation (4.27) 
in |15j . however its derivation is different. As a consequence we expect the stability 
results for intermediate states to confirm earlier known facts. 

4 Stability analysis via the reduced equations 

Determining the (linear) stability of the steady states as given above is an extremely 
simple task in the reduced setting. Indeed, let (9o,7p ) ^ e an equilibrium. The 




(3.3) 



f(J 2 ,cos9) = 0, 
d 2 f(J 2 ,cos9) = 0. 
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linearized equations of motion at this relative equilibrium read as 



" 1 M - -575- (W - Bo) - nmg{R -e cos 9 Q ) 2 9 



T W p(9 )(p 



dtpde 

_& 2 t l 

dtpdO 



(0q)9 — ^mg sin 6o(p, 



where we used that R' n equals mg (|2.6I) at the relative equilibria. It is not hard to 
show that the characteristic polynomial of this system is 



A 3 + [img 



(R — ecos#o) 



T, 



66{VQ 



sin 2 6> 
Tss(^o) 



+ 



(»mg(R ~ e cos Bp) sin ) 2 + (0 g (fl )) ' ^ ( OJ 



V 



'o)TW(^o) 



■o) 



(^mgsin 2 6» O )^(0 O ) 



Tee(0o)?w(^o) 

Due to the translational symmetry in Tp, one eigenvalue is zero. In Appendix [B] 
we show that all remaining eigenvalues have a negative real part, if and only if 
d 2 W/d9 2 {6 ) > 0, or if 



mge cos 9 > W e +c'flcose ( 



/- (A/C- l)sin 2 6» 



+ g cos 6l - 4^ , , , „, . a f in , 3 , (4.1) 

C u (A/C) sm J 8 + (cos 8 — e/R) 2 /' V ' 

with B given by B := (A — C) cos(9q)+C-^. We will further manipulate this equation 
to retrieve the stability results, compare also with |20j . We will retrieve six groups 
depending on how the inertia ratio A/C relates to the eccentricity e/R. Since in 
the literature results have been expressed in terms of the spin of an initial condition 
at a vertical state, we introduce no :— c ^_^ , which is the value of the spin n at 

9 = for a given Jellet J. Similarly, := — c ^ +e j is the spin of the solution with 

Jellet J at 9 — n. Note that for a fixed J these spins are related by uq = —n^^^. 



Vertical spinning state: 
relation (|4.ip yields 



= 0. For the vertical spinning state 9 = 0, the 



--(1--) 

c v r' 



< 



mge 
~C~ 



1 

R 



(4.2) 



It follows that in Group I, the vertical state 9 = is always stable, while for Group 
II and Group III stability requires that 



Kl<n i: =y3=r ( i--) 



(4.3) 



Vertical spinning state: 9 
relation (|4.1I) yields 



(1 



7r. For the vertical spinning state 9 = w, the 



,4 
C 



> 



mge 
~C~ 



R 



(4.4) 
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This condition is never satisfied for Group III, so 8 — ir is unstable; in the case of 



Group I and II, when 4 < (1 + 4), stability requires 



K\>n 2 =:^g^(l + ^. (4.5) 
Note that for tippe tops of Group I and II n\ > n 2 , with n* := 2 v ° r ? ge , The 



equality holds when 4 = | (l + 



C 2 V A 1 i?J • 

Intermediate states. Recall that intermediate states are determined by (13.2[) : 
/(J 2 , cos #o) = 0. Using this condition, the requirement (I4.1[) becomes < g(cos#o) 
where we set 

g (cosfl ):= (t-l)+ $t*r:iZ?S - ^ 

We now prove that g(cos8o) is strictly increasing and changes sign at a bifur- 
cation point for intermediate states. As we already mentioned, a bifurcation point 
along the curve in the (J 2 , cos#)-plane of intermediate states is determined by the 
conditions 

/(J 2 , cos (?)= 0, 
d 2 f(J 2 ,cos8) = + 

4((£ - 1) costf + - cos 2 8) + (cos8- -|) 2 ) = 0. 

An elementary substitution of the first equation into the second shows that the 
function 82/ {J 2 , cos 8) along the intermediate states can also be written as 

J2 

82/ (J 2 , cos 8) = g(cos 



mgeR 2 C 



Hence, bifurcation points are given by those 8 such that <?(cos6V)) = 0. Solving for 
cos 8 gives the two solutions: 

e/R 1 VS 



Note that 1 — A/C — (e/R) 2 > only for tippe tops in Groups I and II. Moreover, 
the solution cos 6b = x-a/c + • • • l ea ds to a contradiction since, for tippe tops of 
Group I, it is incident with the interval ]0, 9 C [ where no intermediate states exist 
and for tippe tops in Group II satisfying 1 — A/C — (e/R) 2 > 0, the number 1 J^y c 
is greater than 1, implying that the + solution of (|4.7I) can not equal a cosine. We 
denote the — solution by Xb, i-e. 



*6 : = T^/c - T=Ajcf VA/C^l - A/C - (e/Rf. 

We conclude that a bifurcation point for intermediate states exists if 1 — A/C — 
(e/R) 2 > and |ar 6 | < 1. 

Before studying in further detail these two conditions, we first show that the 
function g(cos8) is strictly increasing for increasing 8. This result implies that, if 
a bifurcation exists (i.e. a point 8 with g(cos8) = 0) then stability will change. On 
the other hand, if no bifurcation occurs, the entire branch of intermediate states is 
either stable or unstable. 

Let us assume that x = cos 8. If we consider &2f(J 2 ,x) as a function on the 
submanifold f(J 2 ,x) — and if we compute its derivative w.r.t x, i.e. 82.2! — 
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(9\,2f)(d2f /dif), then after some tedious computations we may conclude that the 
sign of this derivative is opposite to the sign of 

8((A/C l)x + e/R) 2 + (A/C l? {A/ ^^ff ? > 

Hence, if there is a bifurcation at a certain x = cos(9b) in the set of intermediate 
states, then we know that the intermediate states for which 9 > 9b are stable, while 
the other branch is unstable. 

It now remains to study the conditions for the bifurcation point to exist. The 
first condition says that 1 — A/C — (e/R) 2 > 0, implying that we only have to 
consider Groups I and II. We start with Group I. 

Group I From (A/C - 1) < -e/R, it follows that 1 - A/C - (e/R) 2 > and 
Xb> < 1. We have to distinguish between two subgroups: Xb < — 1 (Group la) 
and Xb > — 1 (Group lb). From the previous paragraph it should be clear that if 
Xb < — 1 then the value of the function g(x) on the intermediate states — 1 < x < 1 
is negative, i.e. the entire branch is unstable. 



Group II We define three subgroups: Group Ha is the group for which 1 — A/C — 
(e/R) 2 >0and|a;b| < 1, Group lie is defined by 1 - A/C- (e/R) 2 > and x b < -1 
and thirdly Group lib as the group containing the remaining tippe tops in II. Again 
from the previous, we immediately conclude that the entire branch of intermediate 
states is unstable in Group lie. Group lib can alternatively be defined as the group 
containing all tippe tops for which the branch of intermediate state is entirely stable. 
To show this, we first remark that the remaining tippe tops in II are characterized 
by \-A/C-(e/R) 2 > and x b > 1 or I- A/C- (e/R) 2 < 0. If 1- A/C- (e/R) 2 > 
and Xb > 1 then the intermediate branch is entirely stable. If 1 — A/C — (e/R) 2 = 
then Xb > 1 and the intermediate branch is stable. The remaining tippe tops we 
have to consider are characterized by the condition 1 — A/C — (e/R) 2 < 0. To show 
that g(x) > 0, we consider two subcases: (i) if 1 — A/C < then from (|4.6|) it is clear 
that g(x) > 0, (ii) if (e/R) 2 > 1- A/C > then it suffices to compute g(0) (g(x) will 
not change sign since there is no bifurcation point Xb)- from 1 — A/C — (e/R) 2 < 0, 
A/C < 1 and (e/R) 2 < 1 we find 

.9(0) = (A/C 1) + 4 ^g^p > (e/R) 2 > 0. 

The above argument also proves stability for intermediate states in Group III. 

Note that tippe tops in Group II are real 'tippe tops' since they admit tipping 
from a position near 9 = to the inverted state near 9 = n. Tipping never occur for 
tops of Group III, though they may rise up to a (stable) intermediate state. Tippe 
tops of group I never flip over since the position 9 = is always stable. 



4.1 Tippe Top Classification 

The following schematic classification summarizes the previous analysis, and pre- 
sented in Fig. [TJ Compare also with Fig. 3 in [IS] . 

Group I: A/C - 1 < -e/R 

- The non-inverted vertical position 9 — is stable for any value of J. 

- The inverted vertical position 9 — ir is stable for |n w | > 712, unstable otherwise, 
with ri2 given by (|4.5|) . 

- Intermediate states do not exist for all values of 9, but only for 9 > 9 C = 
arccos((^)/(l - £)). 

Group la: - ^^f^A/C^/1 - A/C - (e/R) 2 < -1 . 
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The entire branch of intermediate states is unstable. 
Group lb: if -1 < - -^^^^A/C^JX - A/C - (e/R) 2 := cos(9 b . 

There is a bifurcation: intermediate state are stable if 9 > 9b and unstable if 
9 < 9 b . 

Group II: -e/R < (A/C - 1) < e/R. 

- If |no| > n\ the equilibria 9 — become unstable, with n\ as in (|4.3|) . 

- If In^l > ri2 the equilibria 9 = it become stable. 

- There are intermediate states for any 9. We distinguish the following three sub- 
groups. 

Group Ha: (A/C- 1) < -(e/R) 2 and 



< 1 



T^XJc ~ T^f^fAjC^l-A/C-(e/Rf 

There is a bifurcation of intermediate states. 
Group lib: either (A/C - 1) > -(e/R) 2 or 

T^7c - T^Tc f - A/C (e/R) 2 > 1 

The branch of intermediate states is entirely stable. 
Group lie: (A/C- 1) < -(e/R) 2 and 

e/R 1 </3 



T^jc ~ T^kcfV^V^ A/C -(e/R) 2 < -1 ■ 
The branch of intermediate states is entirely unstable. 
Group III: (A/C - 1) > e/R. 

- The equilibria with 9 = become unstable for |no| > n\. 

- The equilibria with 9 = tt are always unstable. 

- For these tippe tops intermediate states do not exist for 9 g]9 c ,tt[. Bifurcations 
in intermediate states do not occur and the intermediate states are stable. 

Tippe tops in Group II exhibit 'tipping' behavior: if the initial spin satisfies |no| > 
max(ni, ^2 -§37 JE then tipping is possible from 9 = to 9 — n. 



5 Remarks and Conclusions 

We would like to remark that the results on the linear stability for the tippe top 
presented here are equivalent to the stability properties for the model of the tippe 
top without the assumption that the translational friction terms can be neglected. 
In fact, the function g(cos9) which fully determines the stability of the intermediate 
states can be retrieved in the expressions for the eigenvalues of the linearized equa- 
tions of motion about the relative equilibria of the full system. It is also remarkable 
that the presented stability analysis does not depend on the friction coefficient [i, 
which suggests that the above stability analysis is valid for a rather large class of 
possible dissipative friction forces at the point of contact as was pointed out in [15] . 
Up until now, a linear stability analysis for the intermediate states was not available 
upon our knowledge in the literature. 

Most recently Ueda et al. [20 analyzed the motion of the tippe top under the 
gyroscopic balance condition (gbc) and approached the stability problem by per- 
turbing the system around a steady state and obtained under linear approximation 
a first order ode for the perturbation of the variable 9. They derive stability crite- 
ria in terms of the initial spin n given at the non-inverted position 9 = 0. Possible 

3 About the relation between n± and 712: in Group lib n^/nf > 1 holds and in Group lie 
n\/n\ < 1. 
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intermediate states for tippe tops of Group I were not considered. Our approach is 
based on a complete different technique, namely Routhian reduction and we obtain 
a more refined and exhaustive classification of tippe tops in six groups (instead of 
three). Moreover, our analysis does not exclude the possibility of launching the top 
with its stem down (i.e. 9 near n). 
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A Routhian reduction for dissipative systems. 

Assume that a Lagrangian L is given, defined on the tangent space of a manifold 
M on which a local coordinate system (q , . . . ,q n ) is chosen. The system is not 
conservative, in the sense that there is a one form on M, denoted by Q — Qidq 1 
representing the generalized force moments of the non-conservative forces. The 
equations of motion read: 

d fdL\ dL n „ fK . 

dt{w)-d¥ = Qi for i = 1 , . . . , 7i (A.1) 



Theorem 2 (cf. [22]) Assume that there exists a coordinate, say q 1 such that 
(i) J^r = 0, (ii) Qi(q 2 , . . . , q n ) is independent of q 1 for all i = 1, ...,n, (Hi) 
Qi(q 2 , . . . , q n ) = 0; then 8 — dL/dq 1 is a constant of the motion. Assume that the 
latter equation is invertible and allows us to write q 1 = f(q 2 , . . . , q n , q , . . . , q n , 8). 
The system (IA.1I) is equivalent to the system 

d fdR\ OR „ /A oX 

where the new Lagrangian R is defined by R = L — q l 8 and such that all occurrences 
of q 1 in R and Qi are replaced by f. 

By equivalent systems we mean that any solution to (| A. 1 1) with fixed momentum 
dL/dq 1 = 8 is a solution to (|A.2|) and vice versa. 
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B Roots of the characteristic polynomial 

We show that the roots of the polynomial 

f (R - ecos6> ) 2 sin 2 6> 
A + fimg ' 



9o)7V^(6*o) Tgg(9 n ) 

| ( M mg sin 2 60)^(60) 

all have negative real parts if and only if ^g^-(9o) > 0. For our convenience we use 
the following shorthand notations 

^ng(R- e cos 6> ) 2 n a f^ 6 * ) 
a :— — — — — > 0, p ■— 



Tee(0o) ' K ' vWW 

= ^mg sin 2 9 > Q &nd ^ . = W^o) 

Ttp^(9q) Tgg(9o) 

and the polynomial under consideration becomes 

A 3 + {a + 7 )A 2 + (0*7 + /3 2 + <S)A + 57. 

Note that if sin#o = the system (|2.13[) is singular and a coordinate change is 
necessary. The characteristic polynomial written above is not singular. Indeed, we 
have that T^(0o) is proportional to sin 2 9q or that 7 is well-defined in the case 
sin6>o = 0. Similarly, we have that ^J^ g (#0) is proportional to sin#0i and P is well- 
defined. If we now write the necessary and sufficient conditions for this polynomial 
to be Hurwitz |10j . we obtain the conditions 

(a + 7) > 0, (a + 7) (07 + (3 2 + 5) - £7 > and Sj > 0. 

The first condition is trivially satisfied. We now show that the third condition S > 
implies the second condition. We can rewrite the second condition as 

(1 + 077X07 + /3 2 +5) > S, 
and this is valid if S > or if |^r(0o) > 0. 
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